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ABSTRACT 

Stationary density waves rotating at a constant pattern speed Up would produce age gradients across 
spiral arms. We test whether such age gradients are present in M81 by deriving the recent star 
formation histories (SFHs) of 20 regions around one of M81’s grand-design spiral arms. For each 
region, we use resolved stellar populations to determine the SFH by modeling the observed color- 
magnitude diagram (CMD) constructed from archival Hubble Space Telescope (HST) F435W and 
F606W imaging. Although we should be able to detect systematic time delays in our spatially- 
resolved SFHs, we find no evidence of star formation propagation across the spiral arm. Our data 
therefore provide no convincing evidence for a stationary density wave with a single pattern speed 
in M81, and instead favor the scenario of kinematic spiral patterns that are likely driven by tidal 
interactions with the companion galaxies M82 and NGC 3077. 

Keywords: galaxies: spiral - galaxies: star formation - galaxies: ISM - galaxies: individual (M81) - 
galaxies: kinematics and dynamics 

1. INTRODUCTION 

More than half of the galaxies in the Local Uni¬ 
verse show_^piraI_gaHerns highlighted_by young stars 
(e.g.. iNair fc Abrahamll2010t iLintott et al.l f20 1 lh . How¬ 
ever, the origin of these patterns remains an open ques¬ 
tion. Several theories have been suggested to explain 
the underlying mechanisms connecting spiral patterns 
and star formation. Two models in particular are the 
most widely accepte d: the g lobal de nsity wave t heory 
for long-lived arms dLindbladl 119601 : iLin fc Shul 119641 
and the local gravitational instability model for short¬ 
lived transient arms (e.g. , iGoldreich & Lvnden-Belll 
19651 [Julian fc Toomrdll966clSellwood fc Carlberdll984 

Sellwoodll20lH lElmegreenll2011lb The density wave the¬ 
ory explains active star formation in the coherent spiral 
arms as being a result of gas compression by a density 
wave propagating through a galactic disk. In contrast, 
the local gravitational instability model predicts that star 
formation is stochastic, showing no clear age gradients 
across the spiral arm. However, a wide range of observed 
appearances of spiral structures, from grand-design to 
flocculent, implies that a single model may not fully ex¬ 
plain the formation and evolution of all type of spiral 
patterns. For a d etailed review of spiral structures, we 
refer the reader to lDobbs fc Babal (j20l4) and references 
therein. 

Although multiple mechanisms may be responsible for 
the variety of spiral patterns observed, it is straightfor¬ 
ward to test the spiral density wave theory in individual 
systems. Specifically, one should be able to find a system¬ 
atic spatial ordering among SF/gas tracers with different 
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timescales (e.g., Hi for the cold dense gas, 24 pm for ob¬ 
scured stars, and Hct for the young stars) within spiral 
arms that are supported by quasi-static density waves 
( Robertdll969fl . In this model, one assumes that den¬ 
sity waves are propagating through a galactic disk with 
a constant angular pattern speed, and that the corota¬ 
tion radius, R cr , is defined to be where the spiral pattern 
speed becomes the same as the rotational speed of the 
disk. Inside R cr , materials move faster than the den¬ 
sity waves, and vice versa outside R CI . When atomic or 
diffuse gas enters into a spiral arm (i.e., falling into a spi¬ 
ral potential well), it experiences a shock. The resulting 
compression naturally leads to the formation of molecu¬ 
lar clouds, enhancing star formation. Newly formed stars 
disperse the surrounding molecular clouds through stel¬ 
lar feedback and continue to move away from the spiral 
arm. The resulting signature of this process is a spatial 
sequence of cold molecular gas, obscured star formation, 
massive young O/B stars, and evolved stars from the 
upstream to the downstream inside i? cr , and the same 
sequence, but with the opposite angular variation, out¬ 
side R CI . 

With grand-design spiral galaxies, many efforts have 
been made to measure such angular offsets among differ¬ 
ent SF/gas tracers. However, findings and conclusions 
on the angular offsets have conflicted even for th e same 
galaxies (e.g., M51 and M74). While Ibbvle et al.l (12011H 
found no evidence for a sy stematic ordering of d iffer¬ 
ent tracers in these efalaxies. lTamburro et al.l (12008fl and 
lEgusa et al.l (120091) detected the expected systematic an¬ 
gular offsets as a function of galactic radius. However, 
the amplitudes of the measured angular offsets are differ¬ 
ent from each other and thus the estimated timescales for 
the total star formation process es ar e also different. For 
example. iTamburro et all (|2008l) and lEgusa et al.1 ([20091 1 
reported different SF timescales of 1-4 Myr and 5-30 Myr 
based on the angular offset measurements between Hi 
and 24 /jrn and CO and Ha, respectively. A mong vari - 
ous possible causes for the differences. iLouie et all ( l2013t ) 
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Table 1 

A list of Rcr values for M81 found in the literature. Since 
different measurements were made under different assumptions 
about distance to M81, we converted all values to those at the 
distance of 3.8 Mpc. The mean R cr is ~11.24 kpc. 


Reference 

Rcr (kpc) 

Gottesman & Weliachew (1975) 

12.26-13.25 

Rots ('1975') 

12 

Roberts et al. (1975) 

11.5 

_Visser (1980) 

13 

Sakhibov & Smirnov (1987) 

>12.86 

Elmeareen (19891 

9.84 

Lowe et al. (1994) 

10.52 

Westofahl (1998) 

9.8 

Kendall et al.^2008j 

12.67 

Tamburro et al. (2008) 

9.23 

Fene et al. (20141 

9.5 


showed that the discrepancies are mainly caused by dif¬ 
ferences in the choice of SF/gas tracers used to measure 
the angular offsets. They also concluded that CO emis¬ 
sion traces the compressed gas better than Hi 21 cm 
emission, since a portion of HI emission may come from 
the photo-dissociated gas rather than the direct precur¬ 
sor of SF. In addition, 24 /im emission may also be con¬ 
taminated by the underlying older stellar populations, 
which are not associated with the recen t SF activity (e.g., 
iMurphv et al.l[20TTl iLerov et~aLll2012f) . 

In this work, we test the density wave theory by look¬ 
ing for evidence of star formation propagation across the 
spiral arm in M81 using CMD-based SFH analysis, rather 
than using angular offset measurements among different 
SF/gas tracers. This approach avoids any complications 
related to the choice of SF/gas tracers or to the details 
of emission peak finding techniques and angular offset 
measurements from observations of gas emission. 

We carry out this analysis in M81, which is one of the 
largest disk gala xies in the Local Volume. It is nearly 
face-on (i = 59°: lde Blok et afll2008f) and has 2 symmet¬ 
ric grand-design spiral arms with low foreground extinc¬ 
tion (Ay = 0.266: [Schlegell IT998h . Thus, M81 is a good 
laboratory to test theories of spiral formation in detail. 
To compare our SFH analysis with predictions from the 
density wave theory, we adopt a high qual ity HI rotation 
curve from THINGS (Ide Blok et al .1 [2008h and explore a 
variety of pattern speeds reported in the literature for 
M81. The average of the R cr values found in the litera¬ 
ture is ~T1.24 kpc (see Table |T]), and the corresponding 
pattern speed is ~19 km s _1 kpc^ 1 when combined with 
the rotation curve. Throughout this paper, we assume a 
distance of ~3.8 Mpc (i.e., m - M = 27.9), which was 
derived base d on deep resolved ph otometry of the outer 
disk of M81 dWilliams et al.ll2009D . At this distance, an 
angular separation of 1" corresponds to physical distance 
of ~18.5 pc. 

In Section [5] we briefly describe properties of the M81 
imaging data used in this study. Section [3] gives a de¬ 
scription of our methodology for deriving the spatially- 
resolved SFH of spiral arms. In Section [IJ we present 
and discuss our results. We summarize our conclusions 
in Section [5] 


2. DATA 



Figure 1. The HST footprints (green boxes) of 29 fields, which 
cover the entire optical extent of M 81, are superposed on the Hi 
image (THINGS: I Walter et al.1l2008l ). Hi emission is extended far 
beyond optical emission. North is up. East is left. The arrows 
point out the directions toward companions, M82 and NGC 3077. 
Tidal bridges connecting M81 and companions are prominent while 
the tidal features are not seen in the optical. 


We use 29 fields of optical M81 archival data (Pro¬ 
posal ID 10584; Zezas). All fields were imaged in at 
least two filters (F435W and F606W, including some 
F814W observations in outer fields). These fields are 
contiguous and cover the entire optical extent of the 
galaxy. After making bias and flat-field corrections us¬ 
ing the STScI ACS pipeline OPUS, photometry was 
perfo rmed with the ANGST pipeline (Dalcanton et all 
12009 1. The pipeline is optimized for stellar photom¬ 
etry on ACS images using the ACS module with in 
the DOLPHOT photometric package (iDolohinl 12000!) . 
Briefly, the cosmic rays were identified on a combined sin- 
gle drizzled image using the multidrizzle task (PyRAF; 
iKoekemoer et all 12002 ) before measuring flux of indi¬ 
vidual stars. DOLPHOT modifies the Tiny Tim PSF 
(lKristlll995? ) to account for the effects on the PSF shape 
of the telescope temperature changes during orbit. To 
quantify systematic differences between the model and 
true PSF, DOLPHOT determines aperture corrections 
using the most isolated stars in each field. A detailed 
de scription of the photom etry technique can be found 
in iDalcanton et ah (2009). The final catalog utilized in 
this study contains objects flagged as stars with S/N > 
4 in both F435W and F606W filters that pass the sharp¬ 
ness cuts, (sharpp435w + sharpF 606 w ) 2 < 0.075, and 
the crowding cuts, crowdp435w r + crowdp 606 w < 0.6. 
These cuts ensure high quality photometry by selecting 
point-like sources that are not significantly affected by 
crowding dGogarten et al.l 120091) . Using DOLPHOT, we 
also perform artificial star tests to estimate photometric 
completeness and uncertainties as a function of magni¬ 
tude. 
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Fig ure 2. Deprojected image of M81 in Hi (left, from Walter ct al. 2008), 24 fim (right top, from Dale ct al. 2009), and Ha (right bottom, 
from[Hoopes ct al. 200J]). An arrow denotes the direction of galactic rotation when assuming a trailing spiral pattern. On each image, we 
overplot the best-fit logarithmic spiral pattern derived from the HI (red solid line) of the NE arm that we analyze in this work. We also 
mark the reported average corotation radius at 1.24 kpc on the Hi image (dashed circle). 


We use the Hi image from THINGS dWaltcr et al.l 
I2008D to identify the location of M81’s gaseous spiral 
arms. THINGS provides high-quality 21 cm emission 
line maps with an angular resolution of and spectral 
resolution of 2.6 km s -1 . The beam size of the robust 
weighted Hi maps is ~7.5", corresponding to ^140 pc 
at the distance of 3.8 Mpc. Figure Q] shows the Hi map 
overlaid with th e footp ri nts of 29 HS T fields. We also use 
the Her ( KI > NO:ll looncs et al.ll200lD and 24 /im (Spitzer; 
iDale et aLll2009D images to provi de a quick comparison 
between gas/SF tracers in Section l3Jl (see Figure [2]). 

3. THE SPATIALLY RESOLVED STAR FORMATION 
HISTORY 

3.1. Defining Spiral Arms and Selection of Stars 

To define the spiral arms, we first deproject the Hi 
image assumin g a pos ition angle of 330° and inclination 
of 59° (jdc Blok ct al.l[20Q8h (left panel in Figure O and 
then select density maxima positions by picking the cen¬ 
ter of contours within the spiral arms. After translating 
each (x, y) position to polar coordinates (r, 9), we fit 
the selected density maxima positions with a logarithmic 
function of the form ln(r) = ln(a) + b 8. The azimuthal 
angle ( 0 ) increases with the direction of galactic rota¬ 


tion, which is counterclockwise for M81 since we assume 
a trailing spiral pattern. The pattern is almost a perfect 
logarithmic spiral and the derived pitch angle of the HI 
gas arm, i = arctan(|fr|), is ~15°, which is in g ood agree¬ 
ment with previous studies (e.g., iRotsI Il975t iKennicuttl 
1 19811 iLowe et all 11994 iPuerari et al .1120141 )7 We adopt 
this best-fit logarithmic spiral as the peak of the gaseous 
spiral arms. 

M81 has 2 strong spiral arms, separated by ~180°. We 
focus on the NE arm (highlighted with the red solid line 
in Figure [2]), which does not have a strongly streamed 
feature in the HI map compared to the other stretched 
arm (SW arm). This choice should minimize the effects 
of any ongoing tidal interaction on SFH, although the 
SW arm has a similar pitch angle in the Hi gas map. 
The best-fit logarithmic spiral is also overlaid on the 
deprojected 24 /im (right top) and Ha (right bottom) 
images. Although the observed peaks of both Ha and 
24 /im emission are mostly seen downstream of the best- 
fit spiral arm, they are scattered and discrete. This ob¬ 
servational discreteness can cause complexity and ambi¬ 
guity in measuring and interpreting the angular offsets 
among different SF tracers. Note that in both the Hi 
and Ha maps, the spiral pattern is hardly seen in the 
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Figure 3. The 20 analysis regions around the spiral arm are over¬ 
laid on the deprojected Hi map. Individual regions are bounded 
by dashed circles and solid curves. Dashed circles are located at 
~5.6 kpc, ~8.4 kpc, and ~11.2 kpc from the galactic center. To 
clarify our numbering system, we show nominal numbers for re¬ 
gions located at the most leading and most trailing stripes. 


Table 2 

Numbers of stars, 50% photometric completeness limits in F435W 
and F606W bands, the areas covered, and the best-fit A y and 
dAy values for individual regions. 


Region 

N stars 

F435W 

(mag) 

F606W 

(mag) 

Area 

(kpc 2 ) 

Av 

(mag) 

dAv 

(mag) 

i 

14,850 

27.07 

26.93 

2.93 

0.30 

0.9 

2 

22,650 

27.15 

27.02 

3.11 

0.30 

0.8 

3 

18,476 

27.30 

27.26 

3.05 

0.30 

0.4 

4 

6,578 

27.49 

27.48 

3.41 

0.25 

0.1 

5 

16,776 

27.00 

26.84 

3.03 

0.35 

0.5 

6 

15,978 

27.20 

27.07 

3.01 

0.40 

0.8 

7 

16,729 

27.31 

27.28 

2.95 

0.35 

0.5 

8 

6,288 

27.50 

27.49 

3.06 

0.20 

0.3 

9 

17,326 

27.02 

26.90 

3.12 

0.20 

0.8 

10 

12,488 

27.25 

27.15 

2.90 

0.35 

1.1 

11 

12,438 

27.37 

27.34 

2.84 

0.30 

0.6 

12 

4,801 

27.50 

27.51 

2.70 

0.25 

0.2 

13 

11,942 

27.06 

26.89 

3.17 

0.30 

0.9 

14 

13,220 

27.23 

27.12 

2.75 

0.45 

0.8 

15 

7,656 

27.40 

27.35 

2.80 

0.35 

0.4 

16 

3,350 

27.55 

27.57 

2.36 

0.20 

0.3 

17 

12,879 

27.06 

26.90 

3.18 

0.35 

0.6 

18 

14,890 

27.21 

27.13 

2.65 

0.35 

1.0 

19 

6,351 

27.43 

27.39 

2.71 

0.20 

0.6 

20 

2,983 

27.61 

27.58 

2.01 

0.15 

0.6 


inner galactic region (r < 3 kpc), possibly indicating the 
inner Lindblad resonance. In contrast, the 24 /im image 
shows a signature of short spiral arms in the inner galac¬ 
tic region in addition to the grand-design spiral arms, 
as fi rst recognized both in the 3.6 /im and 8 /im images 
bv iKendall ct ah (2008) and explored in more detail by 
iFeng et al.l (12014 ). The driving forces of these short spi¬ 
ral arms is beyond the scope of this paper. 

Figure [3] shows our 20 (5 x 4) analysis regions around 
the best-fit logarithmic spiral arm. Individual regions are 
bounded by dashed circles, which denote galactic radii 
of ~5.6 kpc, ~8.4 kpc, and ~11.2 kpc, and solid curves, 
which define 5 spiral-shaped stripes. These 20 regions 
enable us to obtain the spatially resolved SFHs around 
the spiral arm using sufficient numbers of stars in each 
region, leading to reliable SFHs. We first define 5 spiral¬ 
shaped stripes along the best-fit logarithmic spiral arm 
based on the perpendicular distance from the best-fit spi¬ 
ral arm (from -700 to +700 pc with a constant width of 
280 pc). The total width of 1.4 kpc is wide enough to 
catch the potential age gradients across the spiral arm 
predicted by the density wave theory (see Section PO for 
timescales of interest). Each spiral-shaped stripe is then 
further divided into 4 different bins based on the dis¬ 
tance to the galactic center (from ~2.8 kpc to ~13.8 kpc 
with a constant radial length Ar ~2.8 kpc). The inner¬ 
most and outermost radial boundaries are chosen to limit 
the analysis to the regions where the bulge components 
are negligible and the spiral pattern is clear. By choos¬ 
ing A r ~ 2.8 kpc, the boundary between the third and 
the fourth bins is aligned with the reported average R cr 
of ~11.24 kpc, and thus the fourth bin is likely outside 
R cr . Therefore, in the case of Bp = 19 km s -1 kpc -1 , 
one can avoid the potential ambiguity arising from the 
coexistence of two oppositely propagating SF within a 
single bin. However, for a pattern speed faster (slower) 
than 19 km s -1 kpc -1 , corotation radius lies within the 


third (fourth) radial bin. We assign nominal numbers to 
these 5x4 regions. Smaller numbers are assigned to 
stripes with larger azimuthal angle (i.e., the leading side 
within R cr , but the trailing side beyond R cr ). Within 
each stripe, the numbers increase with increasing galac¬ 
tic radius. 

The numbers of stars, 50% photometric completeness 
limits (derived from artificial star tests), and the areal 
coverage for individual regions are listed in Table [5] All 
regions have almost the same depth and similar areal 
coverage (2-3 kpc 2 ). The observed CMDs of each region 
are shown in Figure [4] Each column represents differ¬ 
ent radial bins, with galactic radius increasing from the 
left to the right columns. Each row represents differ¬ 
ent spiral-shaped stripes. The azimuthal angle decreases 
from top to bottom. For example, when adopting Bp = 
19 km s _1 kpc _1 , in the first 3 columns that are the ra¬ 
dial bins likely inside R cr , top panels are for the leading 
side of the arm and bottom panels are for the trailing 
side of the arm, and vice versa in the last column, which 
is the radial bin likely outside R cr . 


3.2. Deriving the Star Formation Histories 

In all regions, there are sufficient main sequence stars 
for age-dating. The SFH of each region is derived 
using the MATCH package (jDolnhinl l2002f). w hich has 
been applied to man y stud ie s fe.g.. (Williams et, all 
2009t iGoearten et all 120091 iMcQuinn et al.l 20101 


Williams et al. 12011 : Weisz et al. 2011 : M cQ.uinn et al.1 


2012t IWeisz et al] ~ l2014t iLewis et al.l 1201511 7 The fun¬ 

damental principle behind MATCH is to find a set of 
synthetic CMDs that best reproduce the observed CMD, 
using a maximum likelihood method assuming Poisson- 
distributed data. We use the Padova isochrone set 
dGirardi et al.ll2002f) with u pdated models dMarigo et al.1 
l2008t iGirardi et all I2010D for the CMD fitting. In 
the actual fitting process, stars brighter than a 50% 






















































Testing Density Wave Theory in M81 


5 


CO 

o 

co 

Ll 


Increasing galactic radius 


22 

24 

26 

28 

22 

24 

26 • 

28 

22 

24 

26 

28 

22 

24 

26 

28 

22 

24 

26 

28 


Region 1 



Region 2 



Region 3 





m 


Region 4 




Region 5. 


' ,;p/. ■ 




Region 6 





Region 7 




Region 8 




Region 9 


I w- 


Region 10 


Region 1 1 


Region 1 2 



V . v Z'.y.'i: 



N, 





Region 13 


Region 14 


Region 15 


Region 1 6 





,Fw 1 


,■ Y»’-, - 


Region 17 



Region 18 


•* 4' ,/-W> ; 



Region 19 


]■ ..... 


Region 20 


F’W-s&.V- 

O-C-: ...y:.- 
A Pd 


±. ■‘W~. 


-10 1 2-10 1 


-10 1 2-10 1 


F435W - F606W 

Figure 4. Observed CMDs of 20 regions. Each column represents the different radial bins in the spiral arm. The galactic radius increases 
from left to right. Within each radial bin, azimuthal angle decreases from top to bottom panels. All CMDs have significant populations of 
upper main sequence stars, which are critical for constraining recent SFH. 


completeness limit are used. We run MATCH with a 
fixed dist ance of to - M = 27.9, a Kroupa initial mass 
function (iKrounal 120011 1, and a binary fraction of 0.35, 
and then solve for the best-fit foreground extinction 
(Ay), differential extinction (cL4y, internal to the 
galaxy), and linear combination of SFR and [Fe/H] as 
a function of age. We allow the foreground extinction 
to range from 0.05 to 0.55 mag with a step size of 
0.05 mag and the differential extinction to range from 0 
to 1.5 mag with a step size of 0.1 mag. 

We use equally spaced logarithmic time bins between 
log(t/yr) = 6.6 and 10.15 with an increment of 0.15 dex. 


In addition, we require the metallicity to increase mono- 
tonically with time, since the observed CMDs are dom¬ 
inated by the upper main sequence which contains lit¬ 
tle information to constrain metallicity. This metallicity 
constraint enables us to derive a SFH with a physically 
plausible metallicity evolution. We estimate the random 
uncertainties of SFHs by using the Hyb rid Monte Carl o 
(MC) sampling technique described in iDolohinl ( 2013 1. 
By creating random samples following the probability 
density of SFH solutions, the Hybrid MC technique over¬ 
comes the issue seen in uncertainties computed from the 
traditional bootstrap MC resampling, which underesti- 
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Figure 5. Example of the best-fit CMD from MATCH. Upper 
left: the observed CMD. Upper right: the best-fit model CMD. 
Bottom left: the residual of data-model CMD. Brighter colors for 
CMD bins where the model underproduces stars. Bottom right: 
the residual significance (i.e., variance-weighted residual). 


~ i i r _ c -- --i--i_"rr"c 

from ^17-25 km s -1 kpc -1 (Gottesman & Wcliachew 

1975 

: Rots! 19751 Roberts et al,. 19751 iVisser 

198C 

: ISakhibov & Smirnov! 19871 Elmesrreenl 

1989; 

Lowe et all 1994 

Westpfahl 1998 

Kendall et all 

2008; 

Tamburro et all 

2008t Feng et al 

2014|). Given the 


width of our defined spiral arm, we expect stars to 
move through our analysis region in less than ~150 Myr 
at all galactic radii, even where the relative velocity 
between the density wave and the material is modest 
(~30 km s^ 1 ). Therefore, we focus only on the recent 
SFHs. 

Figure [S] shows the SFHs of individual regions during 
the past 150 Myr. In general, the overall SF is weaker at 
the outer two radial bins, most likely due to the low gas 
surface density (i.e., the decline in the gas surface den¬ 
sity with galactic radius). In the outermost radial bin, 
a weaker shock might also be responsible for lower SFR 
around the corotation radius, due to the small relative 
speed between the gas and the density wave, if it exists. 
For clarity, we multiply SFRs of the outermost radial bin 
(last column panels in Figure 0 by 5. The dashed line 
indicates the average SFR over the past 150 Myr in each 
region. All regions show clear evidence of young stars. 


mates the random uncertainties of SFH for time bins with 
low or zero SFR. The systematic uncertainties, which are 
dominated by the uncertainty in stellar evolution models 
(lDolnhinll2012ll . are only important when comparing the 
absolute SFH to other galaxies or simulations, but can 
be safely ignored when comparing regions where CMDs 
have similar depths within a single galaxy, since all re¬ 
gions will be affected similarly. 

Figure [5] presents an example of the best-fit CMD from 
MATCH, along with the observed CMD, the residual af¬ 
ter subtracting the model from the data, and the resid¬ 
ual weighted by the varia nce i n ea ch CM D bin (resid¬ 
ual significance CMD; see iDolnhinl 12002 ). Good agree¬ 
ment between the model CMD and the observed CMD 
indicates that the derived SFH is acceptable. The best 
foreground extinction and differential extinction fits are 
given in the last two columns in Table [2] The mean 
of the best foreground extinction fits is 0.295 mag, and 
it is con s istent with the value of Ay = 0.266 mag from 
iSchleaell ( 1998 1. The mean of the best differential extinc¬ 
tion fits is 0.61 mag. The highest differential extinction 
value (~1 mag) in each stripe is found in the second 
radial bin whe re dust emission i s strongest within the 
arm (Herschel: iBendo ct al.l 1201 Oft . Without considering 
differential extinction, we found the best-fit Ay = 0.45, 
which is cons istent with the value t hat was found for the 
arm region in [Williams ct all (12009! ), in which differential 
extinction was not taken into account. 


4.2. Star Formation Propagation 

Density wave theory predicts ag e gradients across spi¬ 
ral arms. iDobbs fc Pringld (1201011 explored the spatial 
distribution of young star clusters with different ages 
(^2-130 Myr) for four different galaxy models; a galaxy 
with a stationary density wave, a barred galaxy, a floc- 
culent galaxy, and an interacting galaxy. Their model 
galaxy with a stationary density wave shows an obvious 
trend in age of star clusters across its spiral arm, which 
is a clearly differentiated feature from their other model 
galaxies. If M81’s spiral pattern is truly driven by a tra¬ 
ditional density wave, then we expect to find a peak in 
SF that shifts in age as one moves azimuthally away from 
the spiral arm. The direction of the shift in age, and its 
amplitude, depends on the relative angular velocity be¬ 
tween the gas and the density wave, which itself depends 
on galactic radius. 

To explore whether the recent SFHs of M81’s spiral 
arm are consistent with SF propagation predicted by a 
steady density wave with a single pattern speed, we com¬ 
pare the derived SFHs with the expected age trend across 
the spiral arm in each radial bin for a given single pat¬ 
tern speed. The expected timescale (At) for a given 
width of azimuthal angle (Ad) is estimated based on a 
given single pattern speed and the rotation curve ampli¬ 
tude (Vp(r) ss 200- 250 km s” 1 at the radii we analyze; 
Ide Blok et al.lf2008h : 


4. RESULT AND DISCUSSION 
4.1. Recent Star Formation History 

All regions show similar SFHs at ages > 0.2 Gyr 
within the uncertainties, as would be expected if 
old stellar populations have been well mixed az¬ 
imuthally on timescales longer than the dynamical time 
and/or mixed radially t h rough radial m i gration (e.g., 
Sellwood fc Binnevl 120021 lHavwooJ 120081 iRoskar et al'l 
20081) . Fortunately, the timescales in which we are inter¬ 
ested are less than ~0.15 Gyr; M81 has a rotation speed 


A t(r) 


Ad (r) 
fi(r) - Sip 1 


( 1 ) 


where H (r) = YsAfl . 

Figure □ shows the estimated timescale for SF to prop¬ 
agate through the width of each analysis regions as a 
function of galactic radius for 5 different pattern speeds, 
D p = 17, 19, 21, 23, and 25 km s _1 kpc^ 1 . In principle, 
the timescales varies with galactic radius such that it in¬ 
creases toward the R cr , becomes infinity at the R cr , and 
decreases beyond the R cr . The direction of the SF prop- 
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Figure 6. The SFHs over the past 150 Myr. Panel order is the same as Figure[4] For clarity, we multiply SFRs of the outermost radial 
bin (last column panels) by 5. The dashed line marks the average SFR over the past 150 Myr in each region. 


agation outside R cr is expected to be opposite to that in¬ 
side R cr - In the case of 17 p = 19 km s _1 kpc _1 , R cr is well 
aligned with the boundary between our third and fourth 
radial bins, and thus we are able to calculate the average 
timescales to cross our defined spiral arm in all radial 
bins, which are about 32 Myr, 47 Myr, 109 Myr, and 
137 Myr for each radial bin (from inner to outer). These 
timescales are long enough for SF to emerge from dense 
molecular clouds, which is believed to occur on timescales 
of ~5-30 Mvr (e.g ., lElmegreenll2000l : lOstriker et al.ll200lt 
lEeusa et alJl2009f) . However, as we mentioned above, for 
a pattern speed faster (slower) than 19 km s _1 kpc _1 , R cr 
lies within the third (fourth) radial bin. When R cr falls 


within one of our radial bins, we expect two waves of SF 
propagating in opposite senses in that single radial bin. 
Thus, we do not provide the expected timescales for that 
radial bin (third bin for fl p = 21, 23, 25 km s -1 kpc -1 , 
and fourth bin for H p = 17 km s _1 kpc _1 ) due to the 
ambiguity mentioned in Section [3. II 
Because of the uncertainty in the pattern speed and 
the fact that we do not know a priori in which stripe the 
SFR is expected to peak, we need to test for consistency 
with a variety of pattern speeds and possible time shifts 
in the onset of SF. Although our defined 20 regions cover 
the high-density region of the HI gaseous spiral arm, the 
outer edge of our trailing side stripe does not necessarily 
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Figure 7. Predicted timescale of SF propagation across the entire 
5 spiral-shaped stripes as a function of galactic radius, for a vari¬ 
ety of pattern speeds: 17 km s — ^pc -1 (crosses), 19 km s~ 1 kpc - 1 
(filled circles), 21 km s —1 kpc —1 (open squares), 23 km s —x kpc -1 
(pluses), and 25 km s —1 kpc —1 (open diamonds). Shaded boxes de¬ 
note the extents of the four radial bins. For ft p = 19 km s~'kpc 
corotation occurs at ~11.24 kpc, which is well aligned with the 
boundary between the third and fourth radial bins. For ft p = 
17 km s -1 kpc -1 , corotation radius lies within the fourth radial 
bin while corotation radius for the other three pattern speeds (21, 
23, and 25 km s — 1 kpc - 1 ) lies within the third radial bin. 


mark the spot where SF is initiated by a shock. 

In Figure [HJ we therefore compare our derived SFHs 
with the predicted SF propagation for 5 different pattern 
speeds and at least 6 different time shifts in the onset of 
SF to find any systematic trend in the enhanced SF peaks 
across the spiral arm. Each panel presents the SFHs of 5 
adjacent spiral-shaped stripes within a single radial bin, 
with the azimuthal angle decreasing from top to bottom 
(i.e., from leading side to trailing side inside R cr , and vice 
versa outside R cr ); each grayscale panel therefore corre¬ 
sponds to an entire column of Figure [6] Region numbers 
are given in the y-axis for clarity. The SFH of each region 
is normalized to its average SFR over the past 150 Myr to 
look for correlations between SFR enhancements and vi¬ 
able propagating spiral pattern models, and is presented 
in grayscale where darker (lighter) colors indicate higher 
(lower) SFR than the average SFR. 

In each panel, a solid curve traces the ages of the stellar 
populations that are predicted to have formed right at 
the edge of our trailing side spiral-shaped stripe (i.e., a 
zero time shift in the onset of SF) and moved across the 
spiral arm if there were a stationary density wave with a 
given single pattern speed. Dashed (dotted) lines denote 
the positive (negative) time shifts of 10, 30, 50, 70, and 
90 Myr, meaning SF occurred before (after) entering our 
analysis region. 

Comparing the grayscale to the model lines, the en¬ 
hanced SF peaks across the spiral arm show no clear 
systematic trend in all radii to support a steady density 
wave with a single pattern speed. Instead, the recent 
SFHs are fluctuating and stochastic. Our finding is ro¬ 
bust against both adopting a variety of pattern speeds 
(i.e., the different SF propagation timescales found in 
each row of panels) and applying ranges of time shift 
in the onset of SF (i.e., the dashed and dotted lines in 
each panel). For the innermost radial bin, there is a 


weak trend that fits with a time shift of 50 Myr for any 
given pattern speed cases. However, a SFH with much 
higher time resolution would be required to determine 
which pattern speed fits the best with which time shift; 
our SFH time resolution, especially in older age bins, is 
not sufficient enough to constrain this. In addition, the 
amplitudes of some enhanced SF peaks associated with 
the trend are not strong. 

A successful propagation model must also work at all 
radii for a single pattern speed. However, there are no 
clear trends in the outer radial bins that could be ex¬ 
plained by any of these curves with the time shift of 
50 Myr that possibly fits the innermost radial bin. For 
example, in the third radial bin, there is one possible 
trend that shows better agreement either with a zero 
time shift for fl p = 17 km s -1 kpc -1 or with H p = 
19 km s -1 kpc -1 with a time shift of -10 Myr. How¬ 
ever, these trends are too weak to support any specific 
pattern speeds and value of time shift, and are not re¬ 
flected at other radii. Thus, we detect no evidence for 
the age gradients that would have been expected for SF 
propagation by a tradition al spiral density w ave. 

Our results agree with iFovle et all ( 2011 ), who mea¬ 
sured angular offsets between Hi and 24 /jm in 
M81’s spiral arm, and found no systematic ordering 
as a function of galactic radius. Our results also 
agree with other recent studies that found no signifi¬ 
cant observational evidence favoring the density wave 
theory in oth er nearby gra n d-desi gn spiral galaxies. 
For example, iFerreras et al.l (120121 ) found no offsets 
between Ha and NUV-optical color in NGC 4321. 
iMartfnez-Garcfa fe Gonzalez-Lopezliral (120131) also failed 
to detect color gradients in the grand-design two-armed 
spiral galaxies NGC 578, NGC 1703, NGC 4603, NGC 
4939, NGC 6907, and NGC 6951. 

4.3. Tidally induced grand-design spiral arm? 

The absence of a systematic trend in our recent SFHs 
indicates that the grand-design spiral structure of M81 
is not the result of the gas response to long-lived den¬ 
sity waves with a single pattern speed. In fact, sta¬ 
tionary spiral waves have never been reproduced with¬ 
out imposing very specific conditions in numerical sim¬ 
ulations, and spiral pattern s are found to be tran- 
sient/ recurrent features (e.g.. lSellwoodll20TH: iFuiii et al. 

B Wada et~akl 120111: IGrand et al. 2 0121 ; Babim^jiL 
iSellwood fc Carlberd 120141 b iD’Onghia et al. 

recently showed that non-linear growth of self- 
induced spiral structures by swing amplifiication can 
demonstrate ‘apparent’ long-lived spiral arms in their 
simulation of an isolated galaxy. They suggested that 
‘apparent’ long-lived global spiral patterns are merely 
products of connected self-perpetuating local segments, 
which are in local dynamical balance between shear and 
self-gravity. However, the self-perpetuating spiral arms 
that are locally fluctuating with time do not necessarily 
exhibit the systematic radial dependence of the relative 
amplitude of spiral arms that were found in the observed 
mass surface densit y map of M81 (e.g.. lElmegreenlll989l 
iKendall et al.ll2008l ). 

A snapshot of the current bisymmetric spiral patterns 
of M81 might be a result of the superposition of two 
or more m = 2 modes w ith distinct pattern speeds 
(jSellwood & Carlbergll2014D . If M81’s grand-design spi- 
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Figure 8. Comparison between our SFHs and the SF propagation predicted by the steady density wave theory. The comparison was 
carried out in individual radial bins (increasing in galactic radius from left to right) for 5 different pattern speeds (shown in each row): (a) 
Qp = 17 km s _1 kpc _1 , (b) Q p = 19 km s —1 kpc -1 , (c) fl p = 21 km s -1 kpc -1 , (d) = 23 km s — ^-kpc -1 , and (e) S2 P = 25 km s —1 kpc“ 

Within each panel, each row shows the SFH in a given spiral-shaped stripe; the number of the stripe refers to regions defined in Figure [3] 
The SFHs of individual regions are normalized to their average SFR over the past 150 Myr. Grayscale shows the enhanced (darker) or 
depressed (lighter) SF events compared to the average SFR in each region. The solid curve traces the expected dominant stellar ages 
across the spiral arm by assuming the onset of SF at the edge of our trailing side spiral-shaped stripe. The dashed (dotted) curves consider 
positive (negative) time shifts of 10, 30, 50, 70, and 90 Myr. We do not provide the predicted SF propagation in radial bins that are located 
at corotation, where no dominant propagation is expected. In these plots, a propagating spiral density wave would show enhanced SF that 
tracks one of the curves in all four radial bins of a single row corresponding to the true pattern speed. We see no evidence for this expected 
behavior. 
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ral arms consist of the superposition of two or more 
modes lasting 5-10 galactic rotations, then it is rea¬ 
sonable to expect two or more distinct SF propagations 
across the spiral arm. However, we find no such features 
in Figure [8] This does not mean that we can decisively 
rule out this mechanism, since there are possibilities that 
the relevant pattern speeds are out of our search range or 
that the superposition of multiple modes would partially 
erase/overwrite the each other’s age gradients, increasing 
complexity in the resulting SFH. 

As an alternative model for gra nd -desig n 
two-armed galaxies , Too mre fc Toom rc (1972), 
Kormendv fe Normanl (11979H . iHernauistl (1 1 9901) , and 


Bottemal (1200. 'll) argue that all non-barred grand-design 


galaxies with typical rotation curves (i.e., flat rotation 
curves) are likely to be the result of interactions with 
nearby companions since they frequently reside in 
the interacting systems showing tidal features such as 
bridges and tails. These conditions fit M81; it has no 
strong bar, but has companion galaxies. 

Such an interaction is clearly taking place in M81, 
making tidal interactions a compelling mechanism for 
driving its sp iral arms. N umerical simulations (e.g., 
iThomasson fc Donned [l993l ; Yunl 119991 1 suggested tidal 
interactions to exp lain the current HI g as distribution 
in the M81 group (I Appleton ct al.lll98~lf) . According to 
Yun’s simulation (lYunl|1999h , M81 underwent the closest 
passage about 220 Myr and 280 Myr ago by its compan¬ 
ion galaxies M82 and NGC 3077, respec tively . I n addi- 
tion, m any observation al stud ies ('e.g.. iYun. Ho. fc Loj 
1994; |Chandar^t_aL| 120011 iPerez^Gonzalez^t^l 
2006 : Ide Mello et al. 120081 iKonstantopoulos et al 


2009: ISantiago-Cortes. Mavva. &; Rosa-Gonzalez ~ l2010l 


Hoversten et al . 20111 iLim. Hwang, fc Led 120131 ) have 
found evidence for tidally induced SF in both M81 and 
its companion galaxies: ages of compact star clusters and 
star-forming regions range ^400-500 Myr. Although 
we also detected enhanced SF at ages between 200 to 
350 Myr in all 20 regions, the associated systematic 
errors are significant since the stars sensitive to this 
SF peak are near the 50% completeness limit of the 
data, and therefore this result would require additional 
observations to confir m with confidence. __ 

Recent studies (e.g., lOh et al.ll2008l iDobbs et ahll2010l 
iStruck et al.ll201l1 : lOh et al.ll2015fl explored the physical 
properties of simulated galaxies having tidally induced 
symmetric two-armed grand-design spiral arms. These 
tidally excited grand-design spirals are almost logarith¬ 
mic and typically kinematic density waves, which are the 
result of nested elliptical stellar orbits induced by tidal 
perturbation. These kinematic density waves do show 
radial variation of their pattern speed. However, slower 
differential rotation than the material disk enables the 
kinematic spira l s to survive for a ti me scale of ~1 Gyr 
(lOh et al.ll2008l) . iSpeights fc Westpfahll (l2012f) estimated 
the pattern speed of M81’s spiral structure and indeed 
detected its radial dependence, indicating the spiral pat¬ 
tern does actually follow differe ntial rotation. 

Although IDobbs et all (|2010f) modeled M51 with an 
interacting companion (NGC 5195), there are some sim¬ 
ilarities that we can apply to the M81 group as well. 
During detailed comparison of their model galaxy to 
HST observations of M51, they noticed that two pas¬ 
sages by the companion galaxy at different times are 


able to reproduce observed kink/bifurcations in the spi¬ 
ral arm besides a well-defined two-armed spiral pattern. 
A kink/bifurcation is formed where the younger arm, in¬ 
duced by the second passage, and the older arm, induced 
by the first passage, are connected. Since M81 is also 
believed to have experienced two tidal interactions sep¬ 
arated by ^60 Myr with two companion galaxies, M82 
and NGC 3077, one would expect to observe a similar 
feature in M81. In fact, bifurcation is seen in the de- 
projected Hi image (Figure [2J at around the corotation 
radius, but the gas density of the bifurcated arm at outer 
radius is not as high as that of the main NE arm that 
we have analyzed. The SW arm shows more complicated 
features, such as at least two bifurcations. The exis¬ 
tence of the bifurcation indicates the main arm by the 
earlier passage should be prominent by the time of the 
later passage. Therefore, if the bifurcation in M81 is a 
real product of the two separate passages, we are able 
to infer that the spiral structure of M81 ha s survived for 
at least ~280 Myr. IWilliams et all (I2009T ) conducted a 
SFH analysis for their M81-deep (~2 mag deeper than 
our data) field in the outer disk, and also suggested that 
the M81’s spiral arms are at least 100 Myr old. 

5. CONCLUSIONS 

In this paper, we investigated the validity of the sta¬ 
tionary density wave theory to explain the grand-design 
spiral structure in M81. We analyzed resolved stellar 
populations in 20 regions around the logarithmic spiral 
arm, and derived star formation histories of individual 
regions using CMD-fitting. Our approach avoids many 
uncertainties inherent in using discrete SF/gas tracers to 
measure the angular offsets across the spiral arms. 

To test the assumption that the spiral structure of M81 
is driven by a stationary density wave with a single pat¬ 
tern speed, we estimated the timescales for the disk ma¬ 
terial to cross the spiral arm as a function of galactic 
radius based on the relative velocity between the rota¬ 
tion curve and a given single pattern speed ranging from 
17-25 km s _1 kpc _1 . We then compared the predicted SF 
propagation with the potential age gradient imprinted in 
our derived SFHs across the spiral arm. The resulting 
SFHs shows no systematic age gradient across the spiral 
arm at all radii. Rather the SFHs are stochastic over the 
timescales of interest. This result provides convincing 
evidence that the grand-design spiral structure of M81 
is not supported by the stationary density wave with a 
single pattern speed, but instead is likely supported by a 
tidally induced kinematic waves due to the interactions 
with companion galaxies. In addition, we also discuss 
the bifurcation as additional evidence supporting the 
tidally induced kinematic density waves in M81, which 
will wind up slowly and eventually decay. In conclu¬ 
sion, our recent SFH results agree with other studies that 
have found observational evidence against the traditional 
density wave the ory in other nearby grand-design spi¬ 
ral ga laxies (e.g.jFovlc ct all 1201 It iFcrrerasi et all 120121 
iMartmez-Garcfa &: Gonzalez-L6pezlirall2013f) . 
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